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Abstract. In prior work, the authors developed a method of degenerate perturbation 
theory about the Ising limit to derive an effective Hamiltonian describing quantum 
fluctuations in a half-polarized magnetization plateau on the pyrochlore lattice. Here, 
we extend this formulation to an arbitrary lattice of corner sharing simplexes of q 
sites, at a fraction (q — 2k) /q of the saturation magnetization, with < k < q. We 
present explicit effective Hamiltonians for the examples of the checkerboard, kagome, 
and pyrochlore lattices. The consequent ground states in these cases for k = 1 are also 
discussed. 



1. Introduction 

Despite decades of theoretical and experimental work, frustrated quantum magnets 
continue to be an exciting subject for current research both experimentally and 
theoretically. Analytic approaches to such spin systems have, however, not often 
directly confronted experiment. Theoretically, the semi-classical spin-wave theory (i.e 
l/s expansion) has probably been most successful. For instance, quantum corrections 
to the classical staggered magnetization are known to be small even for s = 1/2 in 
the unfrustrated square lattice. The utility of the l/s expansion, however, diminishes 
rapidly as frustration is increased. 

This situation is at its extreme in the class of maximally geometrically frustrated 
structures, which includes the two-dimensional kagome and checkerboard lattices, and 
the three-dimensional pyrochlore lattice. These structures have the common feature 
that they can be decomposed into distinct "corner sharing" simplexes, clusters of spins 
in which all pairs are connected by nearest neighbor bonds, such that the entire lattice 
is covered by these simplexes, and different simplexes share at most one site and no 
bonds. The nearest-neighbor antiferromagnet on this lattice has the property that its 
Hamiltonian can be written entirely in terms of the sum of spins on each simplex. At 
the classical level, this implies a large degeneracy, since any change in configuration 
which keeps this sum constant on each simplex cannot change the energy. This leads 
to considerable technical difficulties in the l/s expansion, which have been addressed 
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in tour-de-force work by Henley and collaborators [1]. An unfortunate outcome of this 
work is that the leading order (0(1/ s )) corrections (spin-wave zero point energy) do 
not fully split the classical degeneracy in most interesting cases. Given the challenges 
already present at 0(1/ s°), it is not surprising that higher order corrections in 1/s, 
which would be necessary to resolve the degeneracy fully and determine the quantum 
ground state for large s, are so far not available. 

In recent work, we have demonstrated that the resolution of the classical degeneracy 
can be understood in an alternate approach, based on perturbation theory about the 
limit of strong Ising exchange anisotropy. That is, we consider the XXZ model 



and carry out perturbation theory in a for the degenerate ground state manifold. This 
is expected to be a reasonable approximation in many problems of interest. First, the 
semiclassical approach demonstrates rather generally that in these models quantum 
fluctuations favor collinear ordered states, in which all spin expectation values are 
aligned along a particular, e.g. S z axis. Choosing Ising anisotropy only selects this 
axis, but does not prejudice the ordering beyond this choice. Second, some of the most 
interesting applications are to magnetization plateaus, which often appear in frustrated 
magnets (e.g in pyrochlores [4]). General arguments (see e.g. [2]) imply that ordering is 
again collinear on such plateaus. Moreover, since such a plateau occurs in a substantial 
applied field, the component of each spin parallel to this field is clearly on average 
larger than its transverse ones. Furthermore, the introduction of Ising anisotropy does 
not modify the symmetry of the Hamiltonian in this case. In [3], we have demonstrated 
explicitly how to carry out easy-axis degenerate perturbation theory for a magnetization 
plateau on the pyrochlore lattice at half the saturation polarization. The result was 
an effective quantum "dimer" Hamiltonian describing the splitting of the degenerate 
manifold of plateau states. 

In this paper, we describe the generalization of these results to the lattices described 
above, for various zero field and plateau states. Our results are obtained under the 
assumption that each corner-sharing simplex contains q sites, and that the Ising ground 
state manifold is comprised of states with k "minority" (S z = —s) and q — k "majority" 
(S z = +s) spins per simplex. The cases mentioned above correspond to q — 3, k — 1 
(kagome lattice at magnetization M = 1/3M S ), q = 4, k = 2 (checkerboard and 
pyrochlore lattices at M — 0), q — 4, k — 1 (checkerboard and pyrochlore lattices 
at M — 1/2M S ). For the k = 1 cases, the corresponding effective Hamiltonians are 
generalized "quantum dimer models" , and we discuss their ground states based on known 
results and simple arguments. The solution for the ground states of the k = 2 models 
(corresponding to zero magnetization) on the pyrochlore and checkerboard lattices are 
left for future work. 

The remainder of the paper is organized as follows. In section 2, we describe the 
derivation of the effective Hamiltonians. In section 3, we give explicit forms for the 
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above-mentioned lattices. In section 4 we discuss the ground states of the effective 
models in the k — 1 cases. 

2. Resume of Method 

In this section, we review the degenerate perturbation theory (DPT) methods developed 
in [3]. Because the fine details are given already in [3], we will highlight only the main 
points and those modifications needed for the more general applications in this paper. 

For any (apart from a set of measure zero) fixed field H in the Ising limit a = 0, the 
ground states of equation 1 comprise a massively degenerate manifold of configurations 
of constant magnetization - "plateau states". Specifically, in these configurations (as 
indicated in the introduction) all spins are aligned along the field axis, i.e Sf = suj, with 
o~i = ±1, and every simplex contains k minority (<7j = —1) spins and q — k majority 
(<7j = +1) spins. We seek an effective Hamiltonian in this degenerate subspace. 

We follow the Brillouin-Wigner formulation of DPT. The ground state wavefunction 
satisfies the exact equation 

oo 
n=0 

where the operator Q = E ^ n (1 —V)7ii. Here Tio = 7~t\ a=0 , and Tli = 7i — Hq. 
Because the resolvent contains the exact energy E, equation (2) is actually a non-linear 
eigenvalue problem. However, to any given order of DPT, E may be expanded in a 
series in a to obtain an equation with a true Hamiltonian form within the degenerate 
manifold. Each factor of Q is at least of 0(a) due to the explicit factor in H±, with 
higher order corrections coming from the expansion of E. Following the strategy of [3], 
we will obtain the lowest order non-constant diagonal and off-diagonal terms in Ti e Q. 
Because we seek only the lowest order term of each type, it is admissible to replace E 
by Eq in Q, a considerable simplification. 

With this replacement, we need only calculate the successive applications of Q 
and 7i\ in equation 2 starting from some arbitary initial state. This is facilitated by 
the simple geometry of the lattice, and by symmetry. The Hamiltonian in equation 1 
has a global U(l) symmetry, corresponding to rotations about the z-axis. As such, 
the z-component of the total magnetization ^ . S* is a conserved quantity at every 
stage of DPT. This leads to some important properties of the resolvent operator 
U= (Ho-Eo)- 1 : 

K- 1 = ^ r y 5f 5/ - H 5/ - E , (3) 

ij 3 

where T^- is the adjacency matrix, i.e. T^- = 1 if i,j are nearest-neighbors and = 
otherwise. We will need the action of 1Z not on states in the ground state manifold, 
but on an arbitrary virtual state reached in a DPT process. We describe these states 
by integers m,j, with < rrij < 2s, such that = Uj(s — mf). Conservation of 
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magnetization enforces £\ rnjUj = for all virtual states. One may readily confirm the 
additional identity 



which results from the constraints on the spin configurations contained on the two 
simplexes sharing site j. Using these two relations, the inverse resolvent can be rewritten 

as 



Note that the resolvent only involves variables on the sites that are modified in the DPT 
process, since all other rrij = 0. With this observation, we are ready to describe the 
calculations. 

We first consider the off-diagonal term, which must cause transitions from one 
plateau state to another. To maintain the constraint requires flipping spins ±s — > =Fs 
in a closed loop of alternating ±s spins. This is because an open string of spin transfer 
changes the number of minority sites on the simplexes at the two ends of the string: 
one suffers a deficit, and the other suffers a surplus. Since Hi is a sum of spin transfer 
operators on links of the lattice, one may readily deduce the order at which the first 
off-diagonal term appears, without calculating it explicitly. The smallest even-length 
non-trivial loop in the lattice (non-retracing walk residing on more than one simplex) 
is of some even length L (L = 6 for pyrochlore and kagome, L = 4 for Checkerboard) 
and surrounds a plaquette of the lattice. This loop contains L/2 minority sites that 
will be converted to majority sites. Each flipping takes 2s operations of Hi, so that 
the off-diagonal term is of order a sL . For large s, this becomes of very high order and 
entirely negligible. It can however be important for moderate values of s. 

As is clear from this discussion, the lowest-order off-diagonal term generates a single 
process: flipping all spins along one length-L loop. Thus it has the general form 



What remains to be calculated is the coefficient Kl(s) multiplying the plaquette flipping 
operator. In [3], an iterative method was developed to calculate this amplitude by 
considering successive applications of Hi on any initial state. Because all the spin- 
flips occurring in any given process are on a single length-L loop, the calculation is 
insensitive to the global structure of the lattice, and depends only upon the length of 
the plaquette. This implies that, for hexagonal plaquettes, the off-diagonal coeffcient is 
identical to that calculated in [3], even on the kagome lattice! The values for K s (s) for 
s = 1/2, 1, 3/2, 2, 5/2 were calculated in [3]. 

For the square plaquette lattices, the calculation is actually significantly simpler. 
At a given stage in the DPT process, the resolvent can be calculated using the number 
of times every (alternating) link has undergone spin transfer up to that step. For the 
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(5) 



^off diagonal = -K L (s)a Ls jJ2 (I IT • • • ITXU • • • U I + h.c.) . (6) 
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square plaquette, there are only two such links, as opposed to three for the hexagonal 
plaquette. After £ link operations, one finds that the resolvent is given simply by 

nj 1 = je(As-e). (7) 

Applying this formula, the off-diagonal coefficient is determined explicitly for any value 
of s: 

As [(2s)!] 2 

We now turn to the diagonal terms in DPT. Unlike the off-diagonal terms above, 
these occur at a fixed order in a independent of s. To understand at what order this 
occurs, first note that a diagonal DPT process of order n, which starts out from some 
plateau state, and returns to it in the end can have at most n sites modified, since 
each site that is modified must be modified again at least once in order to return to 
its original configuration. From the structure of the resolvent discussed above, we say 
that at n th order, the diagonal effective Hamiltonian is a sum of terms, each of which 
involves at most n Ising spin variables. The diagonal effective Hamiltonian therefore 
has the general structure 

n cS [{a t }] = £ £ £ I II ) /o>-i» ■ • • » <0» (9) 

n G„ ai...o„y(ij)gG„ J 

where G n denotes a "graph" of lines connecting the labels a\. . .a n visualized as points, 
i.e. a set of unordered pairs of these labels. Note that the dependence upon the lattice 
geometry enters only through r^-. This function can be derived from direct evaluation 
of equation 2 using the simplified resolvent, as outlined in [3]. 

Remarkably, for any lattice composed of corner-sharing simplexes, all terms in 
equation 9 with n < L are constants in the ground-state manifold. This was 
demonstrated in [3] for the case of the pyrochlore lattice (where L — 6), but the 
arguments follow more generally. We outline the basic ideas, referring to [3] for details. 
First, one can show that all "contractible" diagrams, in which at least one point in 
G n is connected to less than 2 other points, can be reduced to a term of lower order. 
This follows by explicitly summing over the corresponding index (e.g. a n ), using the 
identity in equation 4 and the constant total magnetization. Second, one shows that 
"non-contractible" diagrams with n < L are constant by a different argument. First, 
one notes that the sites a$ for which n^eGn ^a, aj 7^ must form a compact cluster, 
and for n < L this cluster cannot span the smallest non-trivial loop. Summing over 
all clusters satisfying ri(ij)eG„ ^a iaj 7^ 0, one again obtains constant. This is essentially 
because any permutation of the cluster sites a» gives another cluster, and the sum over 
all such permutations can be reduced using the fact that the allowed spin configurations 
of a simplex are all permutations of one reference configuration (with k minority spins 
on q sites). 

The lowest non-constant term in the diagonal effective Hamiltonian therefore comes 
at order n = L. Indeed, by application of the same arguments described above, only 
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one distinct term at this order is non-constant: the one for which G n consists of a single 
loop connecting all labels, i.e. U (ij)eGL ^a t a J = r aia2 r a2a3 • • • r 0i _ l0l T 0j . 01 . This term 
can be explicitly evaluated, and has non-constant contributions only when (a±, ■ ■ • , ai) 
lie sequentially along one of the length L loops of the physical lattice. Evaluating f Gn 
at the spins of these sites leads to the final expression for the diagonal energy function. 



3. Results 



The diagonal effective Hamiltonian obtained above can in all cases be written as a sum 
over plaquettes 

^diag = y^^pfovir ■ ■ ,cvl) , (10) 
v 

where a-pi ■ ■ ■ a-p^ are the L sites ordered sequentially around the plaquette V. We 
choose the function £ to have cyclic symmetry. For the kagome and pyrochlore lattices, 
the plaquettes are hexagonal, bounded by 6 links. For the Checkerboard lattice, the 
plaquettes are square, bounded by 4 links. 

Remarkably, the manipulations in the previous section imply that, as a function 
£p(<7i, • • • , <tl) is the same for all lattices with the same plaquette size (4,6,8 etc.), 
for any k. The differences arise in the allowed plaquette configurations, which depend 
strongly on q and k. This function is most conveniently described by giving the energies 
for all possible plaquette configurations. For the lattices with hexagonal plaquettes we 
find 

£« = fo = 0, fl = £l = _ £2 = _ £5 = ___, (n) 

Js 4 (512s 2 - 256s + 33) a 6 J(l - 8s) 2 s 4 a 6 



£7 = 



16(3 -8s) 2 (4s- 1) 3 ' * 4 8(4s-l) 5 ' 

6 Js 5 (3072s 5 - 4352s 4 + 2432s 3 - 660s 2 + 88s - 5) a 6 

(3-8s) 2 (4s-l) 5 (24s 2 -22s + 5) ' 6 ~ T ' 

Js 4 (128s 2 - 40s + 3) a 6 4Js 6 (5184s 3 - 3600s 2 + 789s - 56) a 6 



8(1 - 4s) 4 (8s - 3) 3(4s - l) 5 (64s 2 - 32s + 3) 2 (72s 2 - 18s + 1) 
Js 4 a 6 



2(l-4s) 2 (8s-3) ' 

The various plaquettes types are shown in figure 1(a). For the pyrochlore k = 2 plateau, 
all 13 plaquette configurations in figure 1(a) are possible. Not surprisingly, the energies 
for plaquette configurations with up and down spins swapped are identical, since in this 
case (k — 2, q — 4) there are equal numbers of minority and majority sites on each 
simplex. Interestingly, this case can be interpreted as an "ice model", in which the two 
down spins on a tetrahedron indicate the locations of protons on the oxygen atom at its 
center. 

For the kagome k — 1 antiferromagnet at 1/3 magnetization and the pyrochlore 
k — 1 antiferromagnet at 1/2 magnetization (the subject of [3]), there are only 5 allowed 
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(a) hexagonal plaquette 



(b) square plaquette 



Figure 1. (Color online) Hexagonal and square plaquette Ising configurations. The 
minority sites are denoted by solid (red) circles. 



plaquette configurations, those with no two neighboring minority sites. These are the 
plaquettes configurations labeled 0, 1, 3, 4, 5. Their energies are the same expressions as 
in 11. 

For lattices with square plaquettes, we find 

Js 3 a 4 



£o — £q — 0, 



Js 3 (64s 2 - 24s + 3) a 4 



£i = £-i 



2(4s - 
8 Js 4 (8s 2 



l) 2 ' 

-5s + l)a 4 



(12) 



2(4s - l) 3 (8s - 1) ' ° (4s -l) 3 (16s 2 - 14s + 3) ' 

The square plaquettes types are described in figure 1(b). The k = 2 (zero magnetization) 
plateau of the checkerboard lattice allows all 6 plaquette types outlined in figure 1(b). 
For k = 1 (1/2 magnetization), only types 0, 1, 3 are allowed. 

As explained in detail in [3], these energies are actually slightly redundant. To see 
this, first define x a as the fraction of plaquettes in the lattice in configuration a. Then 
the energy per plaquette of any state is given by 



"^diag/"^ — ^ ] ^g£ai 



(13) 



where N is the total number of plaquettes. Now one notices that the fractions obey 
two global constraints: first the fractions must sum to J2 a x a = 1, and second, the total 
fraction of minority sites must be -. Denoting the fraction of minority sites in each 
plaquette configuration by M a , the latter constraint implies ^ a x a M a = -. Because of 
these two constraints, we see that it is possible to shift the energies in equation 13 in such 
a way as to change the energy per plaquette by a constant for all allowed configurations: 
£ a — » Co + ciM a , with Co and c\ arbitrary. By appropriate choice of these constants, one 
can tune two plaquette configuration energies to zero. 



Effective Hamiltonians for some highly frustrated magnets 



8 



4. Lowest energy states for the k = 1 effective Hamiltonians 

The plateau states described above can always be mapped to (sometimes overlapping) 
dimer coverings on lattices with sites at the center of each simplex, and links passing 
through every site. Each site in the dimer lattice therefore has q links emanating from 
it. A minority site in the original lattice corresponds to the presence of a dimer on the 
corresponding link in the dimer lattice. With k minority sites, each site in the dimer 
lattice is covered by k dimers. The diagonal term in each effective Hamiltonian simply 
assigns different energies to the various dimer coverings. For the cases with k — 1, which 
correspond to plateau states with non- vanishing magnetization, these models are of the 
form of generalized Quantum Dimer Models (QDMS)[5], though often with a modified 
and more complicated diagonal term. In this section, we use known results for these 
QDMs and simple optimization arguments to analyze the particular examples for which 
we have derived the effective Hamiltonian in the previous section. 

4-1. Checkerboard lattice 

First we address the simplest effective Hamiltonian - the one obtained for the 
checkerboard lattice. The dimer lattice is the square lattice, and so for k = 1 the 
Hilbert space of dimer coverings is the same one of [5]. 

For k — 1, there are only three plaquette configurations possible, so we shift energies 
so that only the energy of the flippable plaquette (type 3) is non-zero and equal to 
V = £ + £ 3 — 2S\. For s > ~ it is straightforward to show that V < 0. For s = | there 
is a divergence, which indicates that the procedure is invalid in that case. For spin- 
1/2, this occurs because the off-diagonal term appears already at order a 2 , whereas the 
diagonal term appears at order a 4 . Because off-diagonal processes are possible already 
at second order, some of the intermediate projection operators in equation 2 must be 
treated more carefully in that case. However, for s = | the off-diagonal term is in any 
case dominant, so the diagonal term can be neglected. Therefore we take V = in this 
case, and for all s > |, we have V < 0. 

Adding the lowest order off-diagonal term, which flips between the two type 3 
plaquette configurations on a given plaquette, we find the lowest order diagonal and 
off-diagonal terms form a QDM of exactly the same form of [5] 

*qdm = vE ( |D) < n| + |a><a| )-< E ( |n)<a| + |n)<D| ) .a*) 

where t = K 4 (s)a 4s J > 0. 

The model 14 has undergone in-depth numerical scrutiny (see [7]), which shows that 
for V/t < +0.6 ± 0.05 one obtains a columnar phase, in which the dimers preferentially 
sit on staggered columns of parallel bonds. For any value of s, we have V/t < 0. 
Therefore, the ground state is always in the columnar phase. 
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4-1.1. Kagome lattice Next we turn to the kagome k = 1 plateau, in which the 
magnetization is | of the full polarization possible. The dimer lattice in this case is 
the honeycomb lattice, and so for k = 1 the Hilbert space of dimer coverings is that of 
the QDM of [8]. The effective Hamiltonian takes on the form 

Wqdm = E^- ( E( |0><0|+|0><01 ) ■ < 15 ) 

■p -p V / 

with the plaquette energies defined in equation 11, and t = K 6 (s)a 6s J > 0. As before, 
for spin s — ~ the diagonal energies are taken to be 0. In that case, the model reduces 
to a particular case of the model considered in [8], in which they find the phase is a 
"plaquette valence bond solid", in which resonating plaquettes (superpositions of the 
two type 5 configurations) form a V3 x a/3 sublattice of all plaqettes. 

We can tune the honeycomb lattice plaquette energies so that types and 3 
have energy 0. The remaining plaquette energies are V 5 = \ (£ + 2£ 5 — 3£ 3 ), V\ = 
\ (— £0 + 2£i — £ 3 ), V4 = £4 — £ 3 . From equation 11, we find that for s > 1, V5 < and 
Vi t 4 > 0, so that the flippable plaquettes (type 5 in figure 1(a)) are favored energetically. 

Ignoring the off-diagonal term for a moment, the lowest energy dimer covering 
configuration for the entire range of s has the \/3 x a/3 structure of the columnar 
state of [8]. This dimer covering includes only type and type 5 plaquettes, with | 
of the plaquettes in the type configuration, and the remaining | plaquettes in the 
type 5 configuration. From the constraint on the total fraction of the minority sites 
J2 a x aM a = I we can immediately deduce that x 5 M 5 < | since all x a > and all 
M a > 0. With M 5 = i we find that X5 < |, so that \/3 x a/3 state has the maximum 
fraction of type 5 (flippable) plaquettes we can pack on the lattice, already an indication 
that this is a very good energy state. Since the remaining plaquettes are of type with 
energy 0, which is the second best energy possible for any plaquette configuration, this 
is clearly the lowest energy possible for any dimer covering in this model. 

Next we will try to analyze the model 15 with the off-diagonal included. Since the 
diagonal term favors flippable plaquettes, it is not unreasonable to approximate ~ 0. 
Then our model becomes the same as the QDM of [8] 

Wqdm = vEf 3XOMOXC '1)-^ (IOXOHOX. ;|) , (16) 

where V = V5 < 0. Calculating the off-diagonal coefficient t by the methods of 
section 2, we obtain a^ s ~^V/t = -0.246927, -2.24511, -45.155, -1228.23 respectively 
for s = 1,3/2,2,5/2. For all s > 1, even when extrapolating a — > 1, we find that the 
ratio V/t takes on a value that puts it in a columnar valence bond solid phase, according 
to the results of [8]. For s = 1, the diagonal and off-diagonal terms are of the same 
order in a, and the ratio V/t ~ —0.247 is a-independent. However, this puts the system 
described by equation 16 within the quoted error bar for the transition point between 
the plaquette and columnar valence bond solid phases in [8]: V/t = —0.2 ± 0.05. The 
s = 1 case may therefore be more sensitive to the values of than other values of s. 
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4-2. Pyrochlore lattice 

Finally we turn to the case of the pyrochlore lattice. The k = 1 case was discussed 
in [3], but we will give a brief summary here. As discussed in section 3, the plaquette 
energies are off-diagonal term are identical to those on the simpler kagome lattice. Thus 
the effective Hamiltonian is the same as the one given in equation 15, but with dimers 
on links of the diamond lattice. For sufficiently large s, the off-diagonal term can be 
neglected, and one need only find the classical dimer covering minimizing the diagonal 
energy. For large s 1, this classical optimization problem could be analytically 
solved. The resulting "trigonaV state has trigonal symmetry and a seven-fold enlarged 
magnetic unit cell. Numerical investigations estimate that this state obtains for s > 2. 

For spin s = 3/2, the situation is more complex. Ignoring first the off-diagonal term, 
the diagonal term alone does not fully split the degeneracy of the dimer manifold, and the 
ground state subspace grows exponentially with the system length. If an extrapolation 
is made to the isotropic limit a — 1, the off-diagonal term is non- negligible. One may 
approximately include its effects by using equation 16, as discussed above for the kagome 
lattice. In this model, two likely ground states were proposed in references [3], [2]: the 
"R" state, with a 4-fold enlarged simple cubic magnetic unit cell, and a "resonating 
plaquette state" , with more complex structure. 

For s = 1 and s = 1/2, the off-diagonal term is dominant. While the above R and 
resonating plaquette states remain candidates, another possibility is a U(l) spin liquid 
state, which is known to be the ground state of equation 16 for V/t close to but less 
than 1. 
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